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ABSTRACT 

The Compton Burst and Transient Source Experiment (BATSE) detected 
hard x-ray flux from Nova Muscae 1991 (GS/GRS 1124-68) during an interval 
of ~100 days beginning ~130 days after the January 1991 main outburst. The 
light curve during this secondary outburst is roughly symmetric, reaching a max- 
imum around mid-July 1991 at an intensity of ~15% of the peak intensity during 
the main outburst. The hard x-ray spectrum displays a soft-to-hard evolution 
during the rise to maximum; the post-maximum spectral evolution is less well 
determined. We compare our observations with those of the GRANAT/SIGMA 
experiment, which covered the initial outburst well but missed most of the sec- 
ondary outburst. 


1. INTRODUCTION 

Nova Muscae 1991 was discovered independently as an x-ray transient 
(GS/GRS 1124-68) in January 1991 by monitors on Ginga (Makino & the Ginga 
team 1991) and GRANAT (Lund & Brandt 1991). In low-energy x-rays ( E < 20 
keV), the source reached a maximum intensity of ~8 Crab on 15 January (Ki- 
tamoto et al. 1992). Subsequently, it decayed almost exponentially with a time 
constant of ~31 days for ~75 days, at which time a secondary outburst occurred, 
followed by a further exponential decay with a slightly longer time constant (~37 
days). The power density spectrum showed flat low-frequency noise (below ~2 
Hz) and variable-frequency quasi-periodic oscillations at higher frequencies (~3- 
10 Hz) (Tanaka, Makino &: Dotani 1991; Grebenev et al. 1992). The spectrum 
was “ultrasoft” at the lowest energies with an apparently power-law tail ex- 
tending to several hundred keV (Gil’fanov et al. 1991). A narrow emission line 
at 481 ± 22 keV was detected during a 13 hour SIGMA observation on 20-21 
January; this was interpreted by Gil’fanov et al. as a red-shifted positron anni- 
hilation line. The power-law tail diminished in intensity faster than the ultrasoft 
component during the first few months. Nova Muscae was not detected in high- 
energy x-rays in April, but was present during SIGMA observations in May and 
August at ~10% of the January maximum. 

An optical nova was detected within the x-ray error circle (Della Valle, 
Jarvis & West 1991a, 1991b). Subsequent monitoring showed it to decay at a 
rate similar to the low-energy x-rays. A binary orbital period of 10.42 hours was 
found, from which a mass function of 3.1 ±0.4 Mg was determined (McClintock, 
Bailyn & Remillard 1992), malting the source a leading black hole candidate. 

The BATSE instrument on the Compton Gamma Ray Observatory 
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(CGRO) functions as an all-sky monitor by detecting Earth occupations of 
low-energy gamma-ray sources (Harmon et al. 1992). CGRO was launched in 
April 1991, too late to observe the main outburst of Nova Muscae. Nevertheless, 
we used the occupation method to search for subsequent emission and found 
that Nova Muscae was detectable by BATSE during an interval of ~100 days, 
beginning just before the May 1991 detection by SIGMA. 

2. OBSERVATIONS 

The BATSE Earth occupation methodology has been described by Har- 
mon et al. (1992). We obtained spectra in 16 energy channels spanning the 
range 20-2000 keV by measuring differences in detector counting rates as Nova 
Muscae rose above or set below the horizon. We fit the data for each occupation 
in each energy channel with a model consisting of a source term superimposed 
on a time-varying background. Potential interference from other sources was 
handled as follows: A catalog of known bright or potentially bright sources is 
maintained as part of BATSE mission operations. All such sources whose lo- 
cation was within 60° of the axis of any detector which also saw Nova Muscae 
were treated as possible interfering sources. If the occupation time of any of the 
possible interfering sources fell within 10 s of a Nova Muscae occupation step, 
that step was eliminated from further analysis. Otherwise, terms were included 
in the fit for all possible interfering sources with an occupation step within the 
fit region (±110 s around each Nova Muscae occupation step). 

We searched for emission from Nova Muscae during the interval 28 April- 
31 October 1991 (TJD 8374-8560). However, during portions of this time, in- 
cluding the interval prior to 5 May, the source was too far out of the CGRO orbit 
plane to produce usable occupations. Nevertheless, it was clear from inspection 
of the initial results that a significant signal was detected from Nova Muscae 
during much of the time when usable occupations were available. In order to 
produce a light curve, we summed the individual occupation count rates over 
3 day intervals and deconvolved the spectra using conventional forward-folding. 
For this purpose we assumed a power-law input spectrum with fixed number 
index a = —2. 

The upper panel of Figure 1 shows the light curve derived from the BATSE 
data. The temporal evolution is characterized by a slow, approximately linear 
rise beginning around 15 May, leading to a single broad, symmetric maximum 
around 13 July, and a decay below our detection limit around 1 September. For 
comparison, we show in the same figure the SIGMA observations (Gil’fanov et 
al. 1991). It can be seen from the figure that the maximum flux during our 
observations was ~15% of the primary maximum. We note that the secondary 
maximum fell approximately midway between the last two SIGMA observations. 

We investigated the evolution of the spectrum during the secondary out- 
burst by summing the data over longer intervals (but never longer than one 
CGRO viewing period, usually two weeks), and then fitting various models to 
the resulting count spectra. A single power-law model was found to be an ad- 
equate fit in all intervals, although alternative models such as optically thin 
thermal bremsstrahlung and Sunyaev-Titarchuk Comptonization could not in 
general be ruled out. In the lower panel of Figure 1 we show the evolution of 
the power-law index during our observations compared with the SIGMA results. 
During the secondary outburst, the hardness increases as the intensity rises to a 
maxim um ; the trend after the maximum is unclear because of limited statistics. 
Gil’fanov et al. (1991) noted that their data were consistent with a clustering 
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Fig. 1. Intensity and 
spectral evolution of the 
entire Nova Muscae out- 
burst, combining both 
BATSE and SIGMA 
data, the latter from 
Gil’fanov et al. 1991. 
The flux is integrated 
over 35-100 keV. 
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around an average index a ~ —2.5 with no trend throughout the outburst. The 
BATSE data, however, show a significantly harder spectrum (a ~ —2) around 
the secondary maximum and possibly an even harder spectrum just before the 
source disappears (although the latter may be subject to systematic errors). 

Recently, a technique for producing images using the BATSE occultation 
data has been developed (Zhang et al. 1993). We generated a number of images 
of the region around Nova Muscae to verify that the latter was in fact the source 
of any emission detected in our occultation fitting, and that all significant nearby 
sources had been accounted for. Samples of these are shown in Figure 2, from 
which one can see that no other source besides Nova Muscae is visible above 
55 keV. Below that energy the only other significant source is the recurrent 
transient pulsar GX 301-2, which was included in our analysis. 

3. DISCUSSION 

Two main theoretical mechanisms have been proposed to explain the out- 
bursts of black-hole X-ray novae: the disk thermal instability model and the 
mass transfer instability model (Hameury et al. 1990 & refs, therein). Chen, 
Livio & Gehrels (1993) attempted to explain the entire light curve by a combi- 
nation of these ideas: a main disk instability outburst is followed by a second 
outburst (in soft x-rays) due to mass transfer by evaporation of the outer lay- 
ers of the secondary by x-ray heating. A third outburst also results from mass 
transfer, this time due to more prolonged heating and expansion of the sec- 
ondary’s convective layer. The second outburst in Nova Muscae was clearly 
seen in late March in soft x-rays but not in hard x-rays (Kitamoto et al. 1992). 
The hard x-ray secondary maximum which we observed would then correspond 
in the above scenario to the third outburst. However, if both of the secondary 
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Fig. 2. Images of a region of the sky near Nova Muscae around the time of 
the secondary maximum (27 June-7 August). The images were produced using 
the radon transform and maximum-entropy inversion (Zhang et al. 1993). The 
energy intervals are 25—55 keV and 55—110 keV for the left and right panels, 
respectively. 

maxima simply represent episodes of increased mass accretion through the disk, 
then it is necessary to explain their completely different spectral characteristics. 
Future models must consider the distinctly different temporal behaviors of the 
hard power-law and ultrasoft spectral components in black-hole x-ray novae. 
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